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ABSTRACT 

Using test-particle simulations, we investigate the temporal dependence of the two-point velocity 
correlation function for charged particles scattering in a time-independent spatially fluctuating mag- 
netic field derived from a three-dimensional isotropic turbulence power spectrum. Such a correlation 
function allowed us to compute the spatial coefficients of diffusion both parallel and perpendicular to 
the average magnetic field. Our simulations confirm the dependence of the perpendicular diffusion 
coefficient on turbulence energy density and particle energy predicted previously by a model for early- 
time charged particle transport. Using the computed diffusion coefficients, we exploit the particle 
velocity autocorrelation to investigate the time-scale over which the particles "decorrelate" from the 
solution to the unperturbed equation of motion. Decorrelation time-scales are evaluated for parallel 
and perpendicular motions, including the drift of the particles from the local magnetic field line. The 
regimes of strong and weak magnetic turbulence are compared for various values of the ratio of the 
particle gyroradius to the correlation length of the magnetic turbulence. Our simulation parameters 
can be applied to energetic particles in the interplanetary space, cosmic rays at the supernova shocks, 
and cosmic-rays transport in the intergalactic medium. 

Subject headings: Physical Data and Processes: turbulence; ISM: cosmic rays, magnetic fields 



1. INTRODUCTION 

In several physical systems containing a fluctuating 
magnetic field in a turbulent plasma, the diffusion of 
charged particles in the direction perpendicular to the 
average magnetic field has been recognized to be impor- 
tant. First, in the heliospheric environment, solar ener- 
getic particles associated with impulsive solar fiares, or 
compact point-like sources, have been observed by multi- 
ple spacecraft orbiting on the heliospheric ecliptic plane, 
e.g., STEREO A/B, SOHO, widely separated in longi- 
tude; this has been interprete d as evidence f or strong 
cross-field diffusion (jWibbcrcnz fc Cand[200l . Signif- 
icant perpendicular diffusion has been also invoked to 
explain recurrent cosmic-ray variations at very high he- 
liospheric latitudes, possibly connected w ith Corotating 
Interaction Regions (|K6ta fc Jokipii| [ l9 98f). Observations 
of solar energetic particles at high heliographic latitude 
made by Ulysses, have al so been interpreted as evidence 
for c ross-field diffusion (jZhang et al.l 120031 iDalla et al.l 
[200l . 

Second; in non-rclativistic coUisionless shocks, e.g., in- 
terplanetary shocks or supernova remnant (SNR) shocks, 
the accelerat ion rate depends on the shock obliquity 
(|JokipiillI987f ). i.e., the angle between the magnetic field 
at the shock and the direction normal to the shock sur- 
face. At a perpendicular shock, since the magnetic field 
lines are frozen with the plasma flow, the transport per- 
pendicular to the magnetic field lines allows the particles 
to remain near the shock enhancing their acceleration. 
Recent multiwavelength campaigns of SNRs (from radio 
up to 7-rays) has not yet constrained the magnetic field 
obliquity at the SNR shock despite a burgeoning number 
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of observational evidences of the magnetic field amplifi- 
cation. 

A single model unifying the transport in the directions 
parallel and perpendicular to the average magnetic field 
is needed to understand the propagation of energetic par- 
ticles in a broad class of environments, i.e. interplanetary 
space, supernova remnant shocks and interstellar space; 
however such a model is still missing. Perpendicular dif- 
fusion has been also studied recently disentangling two 
different contributions: field line random walk and gra- 
dient/curvature drif t of the particle guid i ng ce nter from 
the local field line (jFraschetti k, Jokipiil l201l|) . In the 
limit of weak turbulence and of low ratio of particle gy- 
roradius to the magnetic turbulence coherence length, 
the usual assumption of particle magnetization, i.e., par- 
ticle following field lines, is found to fail even in a simple 
and idealized turbulence geometry, such as the three- 
dimensional isotropic model. For this particular turbu- 
lence, this result must a fortiori hold for high-energy 
particles, because, due to larger gyroradius, high-energy 
particles decorrelate from the initial field line earlier than 
low-energy particles. 

In the present paper we investigate via first-principles, 
using Monte Carlo numerical simulations, the transport 
of charged test particles in the direction perpendicular 
to the average magnetic field at early times, prior to the 
spatial diffusion phase. The magnetic turbulence induces 
a decorrelation of the particle velocity from its unper- 
turbed evolution. The two-point particle velocity corre- 
lation function, used here indistinguishably from parti- 
cle velocity autocorrelation, would be constant in time in 
the absence of field fluctuations. The decorrelation time, 
defined as the time beyond which the Lorentz force expe- 
rienced by the charged particle in the ambient magnetic 
field is uncorrelated with the force initially acting on it, 
depends physically on the magnetic energy density and 
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on the ratio of the gyroradius to the turbulence correla- 
tion length. For a broad range of the parameters studied 
in this paper, the perpendicular particle velocity decor- 
relates within a few gyroperiods, i.e., much earlier than 
the typical diffusion time-scale. 

In the numerical experiments performed in the present 
paper, we relate the time-dependence of the two-point 
parallel and perpendicular particle velocity correlation 
function to the diffusion coefficient as given by th e 
TGK formalism (|Tavlodll921l : IGreenlll951l : lKubd[l957h . 
We explore the regime of strong turbulence and use a 
broad range for the ratio of particle gyroradius to tur- 
bulence correlation length; our results apply both to 
non-relativistic and relativistic particles. Our simula- 
tions confirm the exponential, or Markovian, decay of 
the pitch-angle correlation. However, we find an unex- 
pectedly fast decay of the perpendicular correlation, oc- 
curring within a few gyroperiods, which is not accounted 
for in any phenomenological transport model. 

This work is organized as follows: in Section [2] we de- 
scribe previous models for the velocity autocorrelation; 
in Section [3] we specify the particular turbulence model 
and the parameters for the test particles used in our 
simulations; in Section [3] wc summarize the results for 
the behaviour of diffusion coefficient at strong turbulence 
and high ratio of gyroradius to the turbulence correlation 
length; we also show that our previous weak-turbulence 
model based on the separation between field line random 
walk and gradient /curvature drift is in agreement with 
our simulations; in Section [5] we discuss our results for 
the pitch-angle correlation function and the perpendic- 
ular velocity autocorrelation; in Section [S] we describe 
some astrophysical applications; in Section [7] we summa- 
rize and conclude with the observational implications of 
our results; in Appendixes A and B some details of the 
numerical code are outlined. 

2. SPATIAL DIFFUSION MODELS BASED ON VELOCITY 
AUTOCORRELATION 

The investigation of the particle motion in times prior 
to the diffusion phase can be conveniently based on the 
particle velocity autocorrelation, named Lagrangian as 
it makes use of the instantaneous particle velocity, well- 
defined at times as short as the early ballistic phase. 
Moreover, the velocity autocorrelation offers the most 
direct approach to the long-range correlation in the so- 
called anomalous transport regime (jBakunin 2004). The 
Lagrangian velocity autocorrelation has not been fre- 
quently applied to particle transport for two reasons: 1) 
experimentally it is simpler to measure the Eulerian ve- 
locity autocorrelation, relating two points in space sepa- 
rated by a fixed coordinate distance, than the Lagrangian 
velocity autocorrelation, relating the velocity of the same 
physical particle at two different times; 2) numerical sim- 
ulations require using a large number of particles to ob- 
tain good statistics. 

The random walk of charged particles in a static, i.e., 
time-independent, turbulent magnetic field is commonly 
represented as a sequence of many stochastic events, in- 
depende nt of one an other, i.e., a compound Poissonian 
process (|FteIIe^[l97lh . In a magnetic turbulence, every 
scattering encounter has a probability of not occurring 
exponentially decreasing in time, since a scattering event 
will certainly occur as the time proceeds. An exponential 



decay implies that every scattering event is independent 
from any other and the transport is endowed with com- 
plete lack of memory (Markovian process) . 

An exponential form for the autocorrelation of parti- 
cle velocity v{t) moving in a turbulent fluid {{v(t)v{0)) = 
UqC"*/"^) was first heuristically used bv iTavlo^ ()1921[ ) in 
a seminal discussion of the diffusion coefficient in the 
absence of magnetic field; here T has the meaning of 
characteristic time-scale beyond which the fluid density 
variations or turbulent fluid motions cause a jump of the 
particle velocity to a value uncorrelated with the previ- 
ous value. If magnetic turbulence is included, more re- 
cent theoretical arguments led to the exponential form of 
the particle velocity correlation function in the direction 
along the average magnetic field, or, equivalently, pitch- 
angle correlation function (n(0)n(t)) = 3{vz{0)vz{t))/vQ 
in a static magnetic field ()Earl|[T97^. The expo nential 
form of (/x(t)/i(0)) has been shown (|FormanlfT977[ ) to cor- 
respond to the closed for m for the pitch-angle coefficient 
diffusion in lJokipul (jl966D for a quasi-isotropic scattering. 

We briefly recall some previous approaches to the par- 
ticle velocity autocorrelation to estimate the diffusion co- 
efficients K. In weak turbulence regime and for particle 
isotropically scattering in all three space directions, the 
pitch-an gle correlation is found to decrease exponentially 
in time ( Formanlll977D . i.e., the scattering is Markovian, 
with scattering time-scale T|| related to the parallel mean 
free path (A||) by A|| = fT|| and to the parallel coefficient 
diffusion by k\\ = (w^/3)r|| . In the presence of strong tur- 
bulence, numerical simulations (jCasse et al.ll2002l ) found 
an empirical scaling of the scattering frequency (1/t||) 
with the turbulence magnetic energy density. This scal- 
ing is based on two assumptions: 1) the motion perpen- 
dicular to the local field line, and the eventual decorrela- 
tion from the field line at high rigidity, is dominated by 
the turbulent scales smaller than the gyroradius; 2) the 
decor relation time t± is smaller t han the scattering time 
r|| . In lFraschetti fc Jokipiil ()2011[ ). the scale separation in 
the former assumption 1) of iCasse et al.l ()2002D is found 
to hold only for slab turbulence: the decorrelation from 
the field line is governed by the scales smaller than the 
correlation length, but muc h larger than the gy r oscale 
(it should be noted also that iFraschetti fc Jokipiil ()2011[ ) 
assume a negligible power in the magnetic turbulence at 
the gyroscal e). In contrast, for a 3D-is otropic magnetic 
fluctuation, IFraschetti fc Jokipiil ()2011l ) found that the 
decorrelation from the local fleld line has equal contribu- 
tions from turbulent scales both larger and smaller than 
the correlation length. Secondly, in a strong turbulence 
the fluctuating magnetic field transversal to the average 
field becomes c omparable or large r than Bq and the as- 
s umption 2) in JCasse et al. (' 2002D might be violated. 

IBieber fc Mat thaeus (199^ proposed a model (here- 
after BAM model) for the perpendicular transport of 
high-energy (a few GV for typical interplanetary con- 
ditions at 1 AU) charged particles based on a specific 
ansatz on the form of the perpendicular particle ve- 
locity autocorrelation, by supposing that particles sim- 
ply move along the field l ines. Numerical simulations 
([Giacalone &: Jokipiil 11999) have shown that the stan- 
dard quasi-linear theorv (jjokipiil I1966D provides a dif- 
fusion coefficient in the perpendicular direction larger 
than its numerical estimate at protons energy between 
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1 MeV and 100 MeV, wherea s the diffusion coefficients 
in classical scattering theory (jForman fc Gleesonlll975D 
or the BAM model are too small. A more recent the- 
ory for perpendicular transport closer to numerical find- 
ings th an the two f orcmentioned models is the NLGC 
model (jMatthaeus ct al..,2003i) . For magnetized plasma 
flows , NLGC relies on an approximate relation (jCorrsinl 
|1959[ ) between the Eulerian and the Lagrangian velocity 
autocorrelations. The Corrsin relation assumes a dif- 
fusive nature of the particle displacement between the 
two points in space where the correlation is computed. 
However, such a diffusive assumption has been numeri- 
cally found to hold on ly for purely hydrodynamic flows 
as in iKraichnanI ()1977D but not for MHD flows, as far as 
we are aware. Moreover, the NLGC, being intrinsically 
diffusive, can not apply to non-diffusive magnetic tur- 
bulences geometries: different turbulence geometries are 
found to have di fferent effects on the partic le diffusion 
(see for example iFraschetti fc Jokipil (|2011[ )). In addi- 
tion, the field line random walk relies on a free parame- 
ter (named "a" in lMatthaeus et all ()2003[ )) which is not 
provided by any other auxiliary model and must be em- 
pirically determined from numerical simulations. Within 
the idealized compound diffusion model, where particles 
are strictly tied to magnetic field lines, but scatter along 
them and trace back along the same field line, the tem- 
poral evolution of the perpendicular velocity autocorre- 
lation and the consequent diffusion coefficient have been 

Uy, con fi: 

fusion (|K6ta fc Jokipiill2000( ). 

The velocity autocorrelation {vi{t)vj{0)) is known to 
be related to the symmetric part of the instantaneous 
diffusion tensor Kij{t) by the TGK formalism: Kij(t) = 

Jq d£_{vi{^)vj{0)) . The heuristic and physically moti- 
vated BAM forms of the perpendicular velocity autocor- 
relations at gyroperiod scale are: R±{t) = {Vx{t)Vx{Q)) oc 



determ ined analytically, con firming the expected subdif- 
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(u^/3)e~*/'^-* sin(57i), where il = eBQ/m^c is the rcla- 
tivistic particle gyrofrequency corresponding to the un- 
perturbed field Bq. From that assumption it follows that 
the perpendicular diffusion coefficients (symmetric and 
anti-symmetric parts) are k± — KB^lT±/[l + {nT±)'^] and 
KA = KB(r2rA)^/[l -I- (rir^i)^], where t± and ta are the 
decorrelation time-scales {t± — ta as assumed in the 
B AM model) . S u ch a form of k± generalizes the result 
of iForman et al.l ()1974l ) down to low rigidities (below 4 
GV in interplanetary medium). A comparison of the 
time-scale t± in the two forementioned works allows to 
relate the turbulent power at zero wavenumber to the 
intrinsic spatial coefficient diffusion for field line ran- 
dom walk. In the BAM model, it turns out also that 
k±/k^^ = (t_l/t||)/[1 -I- (flT_L)^], which collapses to the 
billiard ball scattering picture ( Forman fc Gleesonlll975l ) 
if T = T|| = T±. In this paper we explore various regimes 
of the velocity autocorrelation. The result is contrasted 
with the BAM model within broad ranges of rg/Lc and 
cr^, where Vg is the particle gyroradius, the turbulence 
correlation length and = (SB/Bq)"^ the magnetic tur- 
bulence normalized energy density. 

Note that the most general form of the total perpendic- 
ular coefficient of diffusion (symmetric part) comprises 
both the field line meandering part, typically dominant 



contribution at low-rigidity [hmfl), and the departure 
from the local field line (kq) , due to gradient /curvatur e 
drift: = kd + kmflV\\ ()Fraschetti fc Jokipiill201lD . 
We estimate in this paper the correlation time for the mo- 
tion perpendicular to the direction of the average field as 
T± ~ (2/3)r^/Kj^, by using the numerically determined 
value of K± (for the factor (2/3) see Sect. 15. 2p . 

3. NUMERICAL METHOD 

In a series of numerical experiments, we consider a pop- 
ulation of charged test-particles gyrating in a magnetic 
field described as follows: we assume a three-dimensional 
magnetic field of the form B(x) = Bq + (5B(x), with an 
average component Bq = BqBz and a random component 
(5B = 6'B{x, y, z) having a zero mean (((5B(x)) = 0), and 
has a turbulence correlation length Lc- We assume in the 
inertial range a scale-invariant, or Kolmogorov, power- 
spectrum in the three space-dimensions: G{k) oc k~^~'^, 
where k is the wavenumber magnitude, /3 — 5/3 is the 
one-dimensional power-law Kolmogorov index and the 
additional 2 accounts for the dimensionality of the tur- 
bulence (for more details see appendix A). The assump- 
tion of a static magnetic field is reasonable if the particle 
speed largely exceeds the Alfven wave speed. 

We perform numerical integration of equations of 
m otion of charged particles combining the code used 
in IFraschetti fc Melial (|2008) with the prescr iption for 
the turbulence introduced in lGiacalone fc Jokip ii (1994), 
and widely exploited in the last two decades in var- 
ious astro physics contexts: (iGiacalone fc Jokipiil 119991: 
Qin et al.l l200l iBvkov et all 120081: IFraschetti fc Melial 
2008( 1. We determine the particle trajectory by nu- 



merically integrating the equation of motion using the 
Lorentz force determined at the instantaneous particle 
position (see Appendix A for details). 

We follow particle trajectories in a three-dimensional 
spatially unbounded region, since particles escaped from 
a bounded computational domain would be removed 
from the simulation and could artificially modify the esti- 
mate of the instantaneous diffusion coefficient. Charged 
particles arc evolved in various realizations of the mag- 
netic turbulence (see Appendix A for details). The final 
time of our simulation runs is empirically determined as 
the computational time where the asymptotic value for 
the diffusion coefficient is attained. 

The relevant parameters are the ratio of the parti- 
cle gyroradius to the turbulence correlations length, i.e., 
rg/Lc, and the normalized turbulence energy density cr^. 
Therefore our treatment applies to energetic particles in 
various astrophysical environments: from the interplane- 
tary space, to the supernova remnant shocks, and to the 
intergalactic medium. 

4. DIFFUSION COEFFICIENTS 

We report in this section the coefficients of diffusion 
parallel and perpendicular, both symmetric and anti- 
symmetric, resulting from our numerical simulations (see 
Appendix B for details). In Sect. [5] we will discuss the 
autocorrelation of the particle velocity making use of the 
diffusion coefficients computed in the present section to 
estimate the characteristic time-scale for the pitch-angle 
and the perpendicular scattering. 

Fig. ([T]) shows the parallel (km), perpendicular («_[_) 



and ant i- symmetric term {ka 



~ t\yx , see also 
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Fig. 1. — Simulated parallel, anti-symmetric and perpendicu- 
lar terms of the diffusion tensor, in units of Kg, as a function 
of (7^, for rg/Lc = 0.1. For weak turbulence, the classical trans- 
port theory prediction for the diffusion tensor apply (see also text). 
As expected, strong turbulence isotropize the particle transport 
(kx ~ A large power in the turbulence makes the field 

isotropic so no drift can be defi ned and ka cannot be determined 
(see also lGiacalone et al.l II1999I) ) 
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Fig. 2. — Perpendicular diffusion coefficient, in units of kb, as 
a function of cr , for rg/Lc = 0.2. Simulated values are compared 
with the solid line and the dashed line, representing respectively 
the field line random walk and the departure from the local field 
line, or gradient /cur v ature dr i ft diff usion coefficient, analytically 
found in Traschctti & Jokipiil |l20Tll), here FJll, Eqs. (49, 45), 
and the dotted line, representing the quasi-linear result, explicit in 
Eq. (B5) of .Giacalone fc Jokipii. tl99a) . here GJ99. 

Appendix B and Fig. (jBip . right panel, computed as 
Kij ~ {viAxj)) of the diffusion tensor as a function of 
(T^ at fixed particle energy conveniently scaled in units 
of Bohm diffusion coefficient hb = vrg/i. In the weak 
turbulence limit, the behaviour of the diffusion tensor as 
a function of tr^ ~ cr^ and k\\ ~ cr"^) is predicted 
by the standard quasi-linear theory for non-relativistic 
particle ene rgies (see for exampl e the d erivation in the 
appendix of I Giacalone fc Jokipiil (jl999l )'). For rclativis- 
tic particles, such a power law dependence is unchanged. 
The KA^ known to be related to the gradient /curvature 
drift velocity, is c onstant as predicted by t he classical 
scattering theory (jForman fc Gleesonl Il975| ) : in classi- 
cal scattering theory ka/h\\ = {\\\/rg)/{l + (A||/rg)^) 
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Fig. 3. — Perpendicular diffusion coefficient, in units of Kg, as 
a function of rg/Lc, for cr^ = 1.0. Simulated values are compared 
with the solid line and the dashed line, representing respectively 
the field line random walk and the departure from the local field 
line, or gradient /curvature drift diffusion coefficient, analytically 
found in Fraschetti & Jokipii (2011), here FJll, Eqs. (49, 45), 
and the dot ted line, re presenting the qu asi-linear result, explicit in 
Eq. (B5) of lGiacalolie"fc Jokipia lfT999l '), here GJ99. 

at weak turbulence. Thus, k^/ks = + 
1 for cr^ <^ 1, regardless of the particle 
energy (see also lGiacalone et all ()1999() ). 

In the strong turbulence limit shown in Fig. ([T]) the 
particle transport is isotropizcd (k_l = Ky). This finding 
may be due to the particular choice of the turbulence 
power spectrum: different turbulence power spectrum 
can result in a different strong turbulence behaviour. The 
K±_ and K|| merge to the same value at crj ~ 30; the 
particle energy is not expected to affect much the value 
of a1. Wc notice that the convergence of k±_ and to 
Kb holds only for sufficiently large particle energy: for 
Tg/Lc ^0.1, the mean free path is much greater than the 
particle gyroradius so the Bohm diffusion underestimates 
the diffusion coefficients at strong turbulence. 

Fig. ^ focuses on the cr^-dependence of the kj_ 
for r„/Lr. = 0.2 a n d < 1. The resuh 

in iFraschetti &: Jokipiil (|2011[ ) is compared with the 
qu asi-linear result, reported explicitly in Eq. (B5) 
of iGiacalone fc Jokipiil (|1999f ). The best agreement 
with the simulations is f ound using the approach in 
IFraschetti fc Jokipiil (|2011| ). For the 3D-isotropic turbu- 
lence, the only case considered in this paper, the gradi- 
ent/curvature drift contribution is smaller by one order 
of magnitude than the field line random walk contribu- 
tion. 

Fig. ^ focuses on the dependence of k_l on rg/Lc for 
values of Vg/Lc les s than unity. The field-lin e random 
walk computed in IFraschetti fc Jokipiil (j2011h provides 
the closest analytical result to the numerical simulations. 
Notice the increasingly relevant contribution of the gra- 
dient/curvature drift term. This shows that as the par- 
ticle energy increases the departure from the local field 
line becomes relevant and the assumption that a parti- 
cle follows the field line must fail, even for small rg/L,.- 
Notice that Fig. ((31) has been derived only for the case 
of 3D isotropic turbulence; it is not known how the gra- 
dient/curvature drift contribution to kj^ will depend on 
other turbulence models. 
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Fig. 4. — Simulated pitch-angle correlation function for = 0.1 
and various particle energies as a function of if!. Subpanel com- 
pares our simulations with a purely exponential decay having ry 
as characteristic time-scale. At early-time, simulations (in black) 
agree with a purely exponential decay by making use of both 
values of K|| = (t>^/3)T||, reported in this paper (in red) and in 
IGiacalone fc Jokipnl II1999I ') (in blue). 
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Fig. 5. — Pitch-angle correlation function for cr^ = 1.0 and var- 
ious particle energies as a function of tQ. As in Fig. (|4jl, the 
subpanel focuses onto the early-time dependence. Intermediate 
energy density of turbulence can be quite satisfactorily reproduced 
by an exponential decay with as characteristic time-scale (see 
also caption of Fig. ((JJ). 

5. VELOCITY CORRELATION 

In this section we study the time-dependence of the 
Lagrangian particle velocity autocorrelation {vi(t)vj(fi)) . 
The resulting simulations arc compared with models for 
pitch-angle scattering and perpendicular decorrelation. 
The decay time-scale of {vi{t)vj{0)) predicted in these 
models depends on the particle diffusion coefficients par- 
allel and perpendicular to the average magnetic field, 
that we have estimated in the previous section. 

5.1. Pitch-angle scattering 

In Figs. (|4l5l6p . our simulations show that the pitch- 
angle correlation drops to zero as a function of time, or 
tVl, where O = eBo/mjc is the relativistic particle gy- 
rofrequcncy corresponding to the unperturbed field Bq, 
for different values of and for various particle ener- 
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Fig. 6. — Pitch-angle correlation function for cr^ = 10^ and var- 
ious particle energies as a function of tQ. In the subpanel the 
simulation (in black) exhibits a temporal decay much slower than 
a simple exponential (in red). 

gies: the stronger the turbulence and larger the parti- 
cle energy, the faster the correlation decay even within 
a gyroperiod scale. The subpanels in Figs. (|4|5|6|) com- 
pare the early-time decay from our simulations with the 
isotropic scattering form e"*/"^" (|Earil973 ). The scatter- 
ing time r|| is estimated by using r|| = (3/z;^)k||, where 
K|| is computed in our simulation runs and compared 
with the values from IGiacalone fc Jokipiil (jl999( ) . The 
subpanel inside Fig. ([4]) confirms the exponential decay 
of the pitch-angle correlation <^1. The bumps in the 
pitch-angle correlation at multiple integers of gyroperiod 
tg = 27r/^2 suggest that, in its helicoidal trajectory, the 
particle velocity keeps higher or der harmonic s of its per- 
turbed periodic motion (see also iCasse eFall (l200l) '). 

Fig.(I5]) shows that for Vg/Lc — 0.75 the bump at iSl = 
27r is smeared out because particles at such an energy (1 
GeV in the interplanetary medium) can travel a distance 
as large as experiencing therefore several scatterings 
within one gyroperiod (compare also with the uppermost 
panel in Fig. (HI)). 

The subpanel in Fig.(IS]) clearly shows for the first time 
directly in terms of (/i(0)/i(t)) that, for strong turbu- 
lence, the simulated correlation drops to zero signifi- 
cantly slower than e"*/"^" . Such a non-Markovian be- 
haviour is found to hold independently on the particle 
rigidity. We have shown in Fig.([T]) that Ky (and there- 
fore T||) saturates to a constant value for ct^ » 1 at fixed 
energy; therefore, the deviation from the Markovian scat- 
tering shown in the subpanel in Fig. ^ cannot be sim- 
ply accounted for with an ad hoc modification of the cr^- 
dependence of T||. Within the quasi-linear limit, i.e., for 
small magnetic field fluctuations, iFormanI (|1977f ) showed 
that the weighted sum over al l the charac teristic times 
that expresses Ky in Eq.(18) of lEaril (|l97l is equivalent 
to the time-integral up to infinity of the parallel velocity 
autocorrelation. The lowest order term of the Earl's se- 
ries for K|[ , corresponding to the case of purely isotropic 
scattering, can be obtained as time-integral of a simply 
exponential decay correlation function, with character- 
istic time given by ry . Higher order terms, that corre- 
spond to anisotropics in the scattering, would produce 
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Fig. 7. — Perpendicular velocity correlation function for Tg/Lc = 
0.1 as a function of tO,. The three panels correspond to three 
values of turbulence energy density cr : 0.1, 1, 10. In every panel 
our simulations (in black) are compared with BAM model using 
the diffusio n coefficients estimated bo th in the present paper (in 
red) and in lGiacalone &: Jokipiil 119991) (in blue). 

a long-duration tail in the correlation similar to what 
we find in the subpaii el in Fig. ([6|) (see also Eq.(32) et 
seq. in iFormanl (|1977f )). We infer from Fig.([6|) that en- 
hanced turbulence produces a "memory effect" in the 
pitch-angle scattering comparable to the effect of the 
anisotropy terms in the weak turbulence. However, we 
do not perform a fit of our simulations due to the theoret- 
ical uncertainties underlying the higher-order character- 
istic time-scales in the Earl scries. We notice that a pure 
exponential e"*/"^" underestimates the pitch angle corre- 
lation at early-time and therefore the asymptotic value 
of the diffusion coefficient k\\ . Since this paper focuses on 
test particle simulations we will not develop an analytic 
discussion here. 

5.2. Field line decorrelation 

In this section we present our results for the perpen- 
dicular velocity autocorrelation on the gyroperiod scale. 
Fig-® shows that the simulated autocorrelation devi- 
ates from the BAM model ansatz even for relatively small 
values of magnetic fluctuations and of Vg/Lc, within the 
expected regime of vahdity of the BAM modcQ. The two 
curves corresponding to k±_ computed in this paper and 
to the value in Giacalone & Jokipii (1999i) agree in both 
the middle and lower panels. For weak turbulence (lower 
panel), BAM model predicts a much less effective atten- 
uation: the simulated correlation is completely smeared 
out beyond tVl Gtt, or t 3tg. At intermediate fluctu- 
ation level (mid-panel), the exponential suppression oc- 

^ An addition a l nu merical factor is present in 
IBieber fc Matthaeui II1997I') to agree with the QLT limit: 
Qtx = {2/3)rg / kmfl (analogous to Eq. 14 in their paper). 
However, the rate of decrease of the correlation does not depend 
significantly on the particular numerical factor used. We use 
Tj_ = (2/3)rg/Kx to compare our simulations with the BAM 
model. 




Fig. 8. — Perpendicular velocity correlation function for = 1.0 
as a function of if! corresponding to five different values of particle 
energy (rg/Lc = 0.019,0.061,0.2,0.75,4.9). In every panel our 
simulations (in black) are compared with BAM model using the 
diffusio n coefficients estimated bot h in the present paper (in red) 
and in I Giacalone fc Jokipiil II1999I ) (in blue) ; see also caption of 
Fig.Q. 

curs faster; this is due to the fact that large fluctuations 
increase the statistical value of the particle gyrofrequency 
{ft = eB/mjc), computed as an ensemble-average, thus 
reducing the gyroradius (r^ = v/fl), i.e., the instanta- 
neous radius of curvature of the particle trajectory; thus, 
the minimum in the velocity component perpendicular to 
the average field is reached at an earlier time (see Fig. ([7]) 
mid-panel). We note also that the increase of the gy- 
rofrequency fl at larger magnetic fluctuations shortens 
the period of the oscillations, in absolute time. These 
effects arc enhanced at larger turbulence (upper panel 
in Fig.(l7|)), where the oscillations are smeared out and 
the shape approaches a purely exponential decay with a 
surprisingly fast drop. 

The energy dependence of the R±{t) is depicted in 
Fig.([8l). We find a good agreement between our simu- 
lations and the BAM model only at large particle en- 
ergies {E ~ 10 GcV or Tg/Lc = 4.9 for interplanetary 
medium, uppermost panel), where the circular motion of 
the particle is weakly affected by the field line meander- 
ing, occurring at smaller scale. However, as the rigidity 
is reduced, the particle motion resonates with a larger 
range of turbulence scales; the velocity autocorrelation 
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Fig. 9. — Perpendicular antisymmetric velocity correlation func- 
tion for Tg/Lc = 0.1 as a function of tf7 corresponding cr^ = 0.1, 10. 
In every panel our simulation (in black) is compared with BAM 
model (in red) using the diffusion coefficients estimated in the 
present paper (as illustrated in Appendix B). 
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XQ. 

Fig. 10. — Instantaneous perpendicular diffusion coefficient, in 
Kb units, comparing the standard average square displacement (in 
blue) and the TGK (in red) methods with the time integral of 
the Rj^(t) provided by the BAM model (in cyan). We chose here 
rg/Lc = 4.9 and = 1.0, in the BAM model regime validity. 

cannot be retained. 

For completeness, Fig. ([9]) compares our simulations for 
RA{t) with the BAM model. Here ta is estimated as 
TA ~ T± ~ {2/3)rg/K±. Our simulations confirm the 
BAM ansatz for the velocity autocorrelation with the 
assumption t± = ta only for weak turbulence within the 
gyroperiod scale {t ^ 2tg for cr^ = 0.1, compare Fig.(l7|), 
lower and upper panels). 

Note that the instantaneous diffusion coefficient in the 
direction perpendicular to the average field is the time- 
integral of {v±{t)v±{0)), from zero to infinity. If the de- 
cay is sufficiently fast, as shown in Figs. (|7|8p even for low 
values of Tg/Lc and small fluctuations (cr^ 1), a sub- 
stantial correction to the BAM modulated exponential 



is required. Currently no analytic models can describe 
such a rapid drop in the correlation. 
Fig. (|TO)) compares the instantaneous K±{t) from the 

BAM model (K_L(t) = d^R±iO) with our simula- 
tions. Clearly the BAM model is only applicable for 
long times due to sinusoidal factors, which make R{t) as- 
suming positive and negative values (see Fig. (fTO)) ). We 
point out in this paper that the early-time oscillations 
in R(t) cause the k± un derestimate previously found 
(jCiacalone &: Jokipiill 19991 ). Therefore, we conclude that 
a model for the perpendicular diffusion needs to ascertain 
an adequate description of the early-time propagation. 
On the other hand the drawbacks of a diffusion model 
can be identified through a study of the early-time prop- 
agation similarly to the study presented here. 

6. SOME ASTROPHYSICAL APPLICATIONS 

In the preceding sections the charged particle motion 
has been described prior to the diffusion regime in terms 
of the particle velocity autocorrelation. This treatment 
applies to various astrophysical environments: energetic 
particles {E = 0.1 — 10"* McV) in interplanetary space, 
at distance of 1 AU from the Sun, gyrating in an ordered 
magnetic field -Bq ~ 5 nT, with a turbulence < 1 and 
a turbulence correlation length Lc = 10^^ AU; particles 
at energy E ^ 5.1x 10" eV to £■ = 4.1 x 10^^ eV diffusing 
upstream or downstream of the non-relativistic shock of 
supernova remnants, likely in an ordered field Bq ~ 3yuG 
with an amplified turbulence < 10*, and a turbulence 
correlation length 3 pc; particles at energy E = 

1.7 X lO^-* eV to £■ = 1.4 x 10^^ eV transported in the 
turbulent galactic medium with Bq = 3//G, ct^ < 1 and 
Lc ^ 10 pc. 

As an example of application of our simulations, we 
consider the problem of propagation of energetic particles 
in the interplanetary medium detected by a spacecraft 
measuring the in situ magnetic-field. Consider a solar en- 
ergetic particle, i.e., proton, with a kinetic energy T — 15 
MeV in an approximately static solar wind Alfvenic per- 
turbation, released by the CME shock propagating from 
a gradual event. At the location of particle release from 
the acceleration region, if the turbulence can be described 
by a three-dimensional isotropic power spectrum, with 
Lc = 10-3 AU and Bq = 5 nT (rg/Lc = 0.75), Fig. M 
shows that, due to perpendicular diffusion, the correla- 
tion drops on the gyroperiod scale (t^ = 27r/ri ~ 14 sec). 
Likewise, if a CLE proton with T = 10 GeV is released in 
a turbulent interplanetary space with Lc = 10~^ AU and 
Bq = 5 nT [rg/Lc = 4.9), the decorrelation will drop on 
the scale of minutes [tg ~ 152 sec). Therefore, particle 
decorrelation from turbulent magnetic field lines needs 
to be considered in tracing the energetic particle tra- 
jectories in the interplanetary medium. The anisotropic 
phase of the observed flux from solar particle events can 
be strongly affected. 

The above result has also application to the galactic 
cosmic-rays below the "knee" of the cosmic-ray spec- 
trum, thought to be accelerated at the supernova rem- 
nant shock, and escaping into the turbulent interstel- 
lar medium. At those shocks an efficient magnetic field 
amplification has been inferred up to values largely ex- 
ceeding the Rankinc-Hugoniot jump across the shock 
through various independent methods, e.g., from the 
shape of radio synchrotron spectra of energetic electrons 
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(|Revnolds fc Ellisod fT99l) to X —ray rims in the rem- 
nant interior jBerezhko et aJ I2003D . Various explana- 
tions have been proposed in the hterature for such a large 
amplified magnetic field: instability by no n-resonant 
cosmic-rays streaming upstream of the shock (!Bclll l200"^ 
or vortical turbulent mo tion seeding downstream mag- 
netic field amplification ([Giacalone fc Jokipiill2007|) . In 
our simulations, assuming a large magnetic fluctuations, 
i.e., {5B/Bof - 10'', 3 pc and Bq = l^G, 

an energetic proton with energy E = 2.1 x 10^^ eV 
[rgjLc = 0.75), will decorrelate on timescale of 2tg = 92 
yrs (see for example Fig. [H corresponding to cr^ = 1.0, 
panel with Vg/Lc = 0.75). Our simulations do not in- 
clude synchrotron energy losses, negligible for protons, 
but likely to be relevant within a decorrelation time in 
the case of energetic electrons. 

7. SUMMARY AND CONCLUSION 

In this paper we have performed Monte-Carlo simula- 
tions of test particles gyrating in a static spatially tur- 
bulent magnetic field with a three-dimensional isotropic 
power spectrum. First, we have computed the depen- 
dence of the diffusion coefficients in the directions par- 
allel and perpendicular to the average magnetic field on 
the turbulence energy density. For weak turbulence, the 
prediction of our previous model for perpendicular trans- 
port based on the separation between the field line me- 
andering and the gradient / curvature drift from the local 
field line is found to be in better agreement with our 
numerical results than the standard quasi-linear theory. 
For the particular power spectrum that we have consid- 
ered and at fixed particle energy, the drift contribution 
is about one order of magnitude smaller than the con- 
tribution from field line meandering. For strong turbu- 
lence, the diffusion tensor becomes isotropic. Wc have 
also computed the dependence of the diffusion tensor on 
the particle energy, specifically the ratio rg/Lc- Also in 
this case the field-line meandering is found to be in better 
agreement with our numerical results than the standard 
quasi- linear theory. The gradient /curvature drift term 
turns out to be important with the increase of particle 
energy; a drift-dominated diffusion cannot be ruled out 
in different turbulent power spectra. 

Second, we computed the dependence of the particle 
velocity autocorrelation in the directions parallel and 
perpendicular to the average magnetic field on the par- 



ticle energy and on the turbulence energy density. The 
pitch-angle correlation drops exponentially in time for 
weak fluctuations, as predicted for the isotropic scat- 
tering case. We also found that for the case of strong 
turbulence the decay is not exponential. The deviation 
from the exponential decay in strong turbulence cannot 
be accounted for by current models. The perpendic- 
ular velocity autocorrelation decays faster than an ex- 
ponentially modulated oscillation predicted by previous 
models. Even in the weak turbulence case and for gyro- 
radius smaller than the turbulence correlation scale, no 
significant correlation is found beyond three gyroperiods, 
whereas the modulated exponential decays much slower. 
Although in a strong turbulence the particle gyroradius 
cannot be uniquely defined, our simulations show that 
the statistical effect of strong turbulence is reducing the 
instantaneous radius of curvature of the particle orbit; 
the correlation is lost within a fraction of gyropcriod. 
However, for gyroradii larger than the correlation length 
the exponential modulation agrees with the simulations: 
on this scale the effect of the turbulence is not relevant 
and the particle pursues a quasi-helicoidal motion in a 
uniform field. However, for smaller gyroradii, the esti- 
mate of the decorrelation time and the magnitude shows 
the need of a new model. 

We do not provide here a phenomenological fit or ana- 
lytical model for the velocity autocorrelation. Our work 
shows that the underestimate of the diffusion coefficient 
previously found in the classical scattering theory can 
be explained in terms of lack of consistent model for the 
scattering and decorrelation times, and sheds light for 
the future investigations. A model for the diffusion co- 
efficient based on the velocity autocorrelation needs to 
provide accurate description of the early-time transport, 
because the contributions to the diffusion beyond the gy- 
ropcriod scale average out to zero. 
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ported by NSF grant ATM0447354 and by NASA grants 
NNX07AH19G and NNX10AF24G. FF thanks S. DaUa 
for useful correspondence. FF wishes to acknowledge 
the hospitality of Observatory of Paris-Meudon (France), 
where part of this work was completed. 



APPENDIX 
A. NUMERICAL SET-UP 

The global magnetic field is written as a sum of a background field Bq, constant and statistically uniform, and a 
turbulent field varying in space and independent on time. The equation of motion of a test-particle with charge e and 
mass m moving with velocity v(t) in a magnetic field B(x) is the Lorcntz equation 

^ = u(t) X fl(x) , (Al) 



where r2(x) = eB(x)/(mc7) with 7 = 1/ ^Jl — (w/c)^ the Lorentz factor, c the speed of light in vacuum and where 
u(t) = 7v(t)/c. We calculate the trajectory of the particle in a magnetic field as a solution of Ea. (|Aip where t is 
the time in the rest frame of the propagation region. The quantity f2(x) in Ea. (|Aip is given by r2(x) = fJo + (5Jl(x) 
where fio = (e/TOC7)Bo, in terms of the background magnetic field Bq, and (5Jl(x) is the turbulent magnetic field. We 
verified that the particle energy is conserved with a relative accuracy of 10"^. We ignore any large-scale background 
electric fields. 

The fiuctuating field comprises of equal intensity circularly polarized left-handed and right-handed Alfven waves. 
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Fig. B1. — Left As a function of time, the diagram compares the coefficients of diffusion kx_ and ft|| computed by using ((Aa::)^)/2t (in 

dashed) and /q (fa: (0))d5i oi' TGK (in soUd) for rg/Lc = 0.2 and tr^ = 1.0. If the TGK formula is used, the number of particles 
necessary to reduce numerical fluctuations is A^p ~ 50, 000. Right As a function of time, the diagram compares the anti-diagonal terms of 



the diffusion tensor, i.e., n^y and Kyx, in units of Kg, anti-symmetric as expected (k.^ 



-Kyx)- Here Vg/Lc 



0.1, 



0.316 and the 



number of particles used is Np = 50, 000. We draw the time-averages to clearly show the asymptotic value. 

The procedure of building the turbulence calls for the random generation of a given number N of transverse waves 
kj, i ~ 1,..,N at every point of physical space where the particle i s found, each with a rand om amplitude, phase 
and orientation. We use the form of the power spectrum provided in lGiacalone fc Jokipiil (|1999( ). Eq.s (3) - (7). The 
wavenumber of the A'^ wave modes ranges from kj^^ji = lir j L^^x fo kj^ax — '^'^/-^min-j 

spanning over four decades 

LmaxI L„iin ~ 10'^. The coherence scale Lc, or bend-over scale, for the turbulence power spectrum, is chosen such that 
Lc = 10^ Lmin- The values of fci-magnitude are logarithmically equispaced, Ak/k = constant. We sampled the power 
spectrum over a discrete number of modes Nm- We choose = 188 modes of magnitude fc^; the resulting coefficients 
of diffusion do not depend significantly on the sampling resolution in fc— space for Nm > N^- 

Wc inject a large number of particles at the same point of space, i.e., the origin of the coordinate system, with 
randomly oriented velocity, but fixed energy. We integrate the Eq. (|Aip by using a time-step adjustable Burlisch- 
Stoer method (iPress et al.lll986[) . which in our numerical solution of 2'"^ order ODE resulted more accurate than a 
Runge-Kutta 5*'' order. This approach is well suited to investigate the particle transport in turbulence. A different 
approach of computing different realizations of the turbulent magnetic field in every cell of a structured Cartesian 
grid is much more time-consuming. The latter approach would also require adapting the lattice spacing in order to 
mainta in the same space resolutio n in physical space. 

As in lGiacalone fc Jokipiil (jl999f ) the randomization is performed over the initial particle velocity orientation and over 
the turbulence realization, shuffled every 50 particles. Including the fluctuating field statistics ensures a meaningful 
comparison with a theoretically computed turbulence power spectrum, which is by deflnition an average over an 
ensemble of field realizations. 

B. COMPUTATION OF THE DIFFUSION COEFFICIENTS 

We computed the parallel (k||) and symmetric perpendicular {k±) coefficients of diffusion in the previously specified 
turbulence using two different methods. The first method uses the formula for the instantaneous coefficient diffusion 
based on the ensemble-averaged mean square displacement: 



dUt) 



1 



^ n—l 



{Xjjt) - Xi{to)){Xj{t) - Xjjto)) 
t-to 



(Bl) 



where Np is the total number of particles (the turbulence realization being shuffled every 50 particles), Xi{t) the particle 
trajectory along the i-th space coordinate and to injection time. 
The second method involves the particle velocity autocorrelation given by the TGK formula: 



1 d 



4-{t) = i^j^{^x,^Xj) 



(B2) 



where Vk{t) is particle velocity along the space coordinate Xk- The last equivalence is valid for fluctuating velocity 
statistically homogeneous for sufficiently long time. 

In both cases the coefficient of diffusion is estimated as the apparent asymptotic value in the diagram representing 
our numerical simulations (see Fig. (|B1|) . Left panel), that is attained after a large number of scatterings off the 
magnetic disturbance: dii{t) kh for t — >■ oo. For Vg <^ Lc it takes many gyroperiods before the force acting upon 
the particle dccorrclate from the initial force. For rg > Lc, the force decorrelates much earlier, allowing eventually 
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the particle to gyrate around a field line different from the field line it was released on. Thus the diffusion time-scale 
ranges from lO^tfi (for Vg <C Lc ) to lO'^tQ (for Vg > Lc). 

By u sing the first method fEa lBl[) . we found results numerically convergent and consistent with lGiacalone fc Jokipil 
(|1999[ ) by injecting only Np ~ 2,500 particles. However, the second method (Eq lB2|l . due to the large numerical 
fluctuations around zero of the velocity correlation at t ^ requires a m uch larger number of particles, i.e., 

Np ~ 50,000, for the results of the two methods to be comparable (see Fig. (|Bl[) . Left panel). The chosen number of 
particles compromises between reducing the computation time and smearing out the numerical fluctuations in the tail 
at m>10- 10^ in Figs. pl5l6)) . 

Equations (|B11 IB2p are symmetric in the indexes i — ?', therefore no t suita ble to define the anti-symmetric part of 
the diffusion tensor. We follow here the argument in IGiacalone et aTl (jl999f ). valid for a nearly isotropic population 
of non-rclativistic particles streaming without convection, to write the anti-symmetric part of the diffusion tensor: 
j), where Vi is the velocity along the i direction and Ax^ the space displacement along the j direction. In 
Fig. (|Bip . Right panel, the averages of Hxy and Ky^ — --Kxy are plotted as a function of time, in order to smear out the 
large fluctuation due to velocity oscillations. 
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